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The standard procedure to determine (analytically) the values of the quark masses is to relate 

QCD two-point functions to experimental data in the framework of QCD sum rules. In the case 
^ ' of the light quark sector, the ideal Green function is the pseudoscalar correlator which involves 

'nI" , the quark masses as an overall multiplicative factor. For the past thirty years this method 

has been affected by systematic uncertainties originating in the hadronic resonance sector, thus 
^f) • limiting the accuracy of the results. Recently, a major breakthrough has been made allowing 

ly^ . for a considerable reduction of these systematic uncertainties and leading to light quark masses 

O I accurate to better than 8%. This procedure will be described in this talk for the up-, down-, 

strange-quark masses, after a general introduction to the method of QCD sum rules. 
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1 Introduction 



Due to quark and gluon confinement in QCD the (analytical) determination of the values of 
the quark masses necessarily calls for an approach different from that for ordinary, non confined 
particles. The ideal approach is to consider a QCD correlation function which on the one hand 
involves the quark masses and other QCD parameters, and on the other hand it involves a mea- 
surable (hadronic) spectral function. Given a framework to relate both representations of this 
correlator, i.e. the QCD and the hadronic representation, the quark masses would then become 
a function of QCD parameters, e.g. the strong coupling, some vacuum condensates reflecting 
confinement, etc., and measurable hadronic parameters. This framework has traditionally been 
that of the QCD sum rules [T], to be described in the following sections. Unfortunately, the 
directly measurable Green functions in the light quark sector are the vector and axial-vector 
correlators entering r-decays, and which involve the quark masses as sub-leading contributions. 
This feature impacts negatively on the accuracy of the results. The pseudoscalar correlator, 
which contains the quark masses as an overall multiplicative factor, thus making it the ideal 
object to be used, is not realistically measurable beyond the ground state pseudoscalar meson 
pole (pion or kaon) . In the case of the scalar correlator the situation is worsened by the absence 
of a ground state pole. The standard approach for close to thirty years in connection with 
the (pseudo) scalar correlator has been to model the hadronic resonance spectral functions as 
best as possible. In the early days, radial excitations of the ground state pseudoscalar meson 
were parametrized in zero-width with model dependent couplings, followed later by finite-width 
parameterizations, and eventually incorporating sound threshold constraints from chiral pertur- 
bation theory [2]-[5]. At present, the masses and widths of the first two radial excitations of 
the pion and the kaon are well known experimentally. However, the model dependency is still 
unavoidable as inelasticity, non-resonant background and resonance interference are realistically 
impossible to guess. 

As a result of this model dependency of the hadronic sector, light quark mass determinations 
have been historically affected by systematic uncertainties. A major breakthrough has been 
achieved recently [6]- [7] by using a successful procedure which reduces considerably the contri- 
bution of the unknown resonance sector. This quenching is accompanied by a corresponding 
enhancement of the better known terms in the correlator, i.e. perturbative QCD and the pseu- 
doscalar meson pole. As a result, the light quark masses can now be determined (analytically) 
with an unprecedented accuracy of better than 8%, as will be described in this talk. 



2 QCD Sum Rules 



The method of QCD sum rules, introduced by Shifman, Vainshtein, and Zakharov [8] more than 
thirty years ago, has become a powerful technique to study hadronic physics in the low energy 
resonance region by means of QCD [1]. This method is also a complementary tool to numerical 
simulations of QCD on a lattice [9]. The range of applications has steadily grown over the years 
and covers the determination of the hadronic spectrum (masses, couplings, and widths), as well 
as electromagnetic, weak, and strong form factors. QCD sum rules have also been extended 
to finite temperature and density [10] , thus allowing for a study of chiral symmetry restoration 
and quark-gluon deconfinement [11]. On the QCD sector, sum rules are employed to extract the 



values of the quark masses [T] and of the strong couphng at a scale of the r-lepton mass [12j- 
|14j . and they are also a tool to confront QCD predictions with experimental data, e.g. from 
e^e~ annihilation and hadronic r-lepton decays. This method is based on two fundamental 
pillars: (i) the operator product expansion of current correlators at short distances, extended 
beyond perturbation theory to incorporate quark-gluon confinement, and (ii) Cauchy's theorem 
in the complex energy (squared) plane, often referred to as quark-hadron duality. To be more 
specific, let us consider a typical object in QCD in the form of the two-point function, or current 
correlator 

n(g2) = i I d^x e^"^ < 0| T{J{x) J(0)) |0 >, (1) 

where the local current J{x) is built from the quark and gluon fields entering the QCD La- 
grangian, and it has definite quantum numbers. Alternatively, this current can be written in 
terms of hadronic fields having the same quantum numbers. The specific choice of current will 
depend on the application one has in mind. For instance, if one is interested in determining the 
hadronic properties of the /?"'"-meson, then one would choose the QCD vector isovector current 
•Jij,{x) = d{x)^fj,u{x), and its hadronic realization in terms of the /9+-meson field. If the goal 
is to determine the values of the light quark masses, then the ideal object would be the corre- 
lator involving the axial-vector current divergences J^{x)Ya = {rrii + mj)ip^{x)j5ipj{x), with i,j 
the up, down, or strange quark flavors. The hadronic representation of this correlator contains 
the pseudoscalar meson (vr or K) mass and coupling, followed by its radial excitations and the 
hadronic continuum. The tool to relate these two representations is Cauchy's theorem in the 
complex energy (squared) plane, to be discussed shortly. 

The QCD correlator, Eq.(l), will contain a perturbative piece (PQCD), computed up to a given 
loop order in perturbation theory, and a non perturbative part mostly reflecting quark-gluon 
confinement. The leading order in PQCD is shown in Fig.l. Since QCD has never been solved 
analytically, the effects due to confinement can only be introduced by parameterizing quark and 
gluon propagator corrections effectively in terms of vacuum condensates. This is done as follows. 
In the case of the quark propagator 

SFip) = -^ =^ W^' (2) 

p — m p — m + 2^(p ) 

the quark propagator correction S(p^) would contain the information on confinement. One 
expects this correction to peak at and near the quark mass-shell, i.e. for p ~ in the case of 
light quarks. Effectively, this can be viewed as in Fig. 2, where the (infrared) quarks in the 
loop have zero momentum and interact strongly with the physical QCD vacuum. This effect is 
then parameterized in terms of the quark condensate (0|^(0)g(0)|0). Similarly, in the case of the 
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Figure 1: Leading order PQCD correlator. All values of the four-momentum of the quark in the 
loop are allowed. The blue wiggly line represents the current of momentum q. 



gluon propagator one would have 

where the gluon propagator correction will peak at fc ~ 0, and the effect of confinement in this 
case will be parameterized by the gluon condensate (0|q;s G^'^ ■ G^,/|0) (see Fig. 3). In addition 
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Figure 2: Quark propagator modification due to (infrared) quarks interacting with the physical 
QCD vacuum, and involving the quark condensate. Large momentum fiows only through the 
bottom propagator. 

to the quark and the gluon condensate there is a plethora of higher order condensates entering 
the OPE of the current correlator at short distances, i.e. 

n(g')lQCD = Co/ + ^ C2^+2(g',/u2) (0|d2iV+2(;u')|0) , (4) 

Ar=0 

where [i^ is the renormalization scale, and where the Wilson coefficients in this expansion de- 
pend on the Lorentz indexes and quantum numbers of J{x) and of the local gauge invariant 
operators O^ built from the quark and gluon fields. These operators are ordered by increasing 
dimensionality and the Wilson coefficients, calculable in PQCD, fall off by corresponding powers 
of —q^. Since there are no gauge invariant operators of dimension d = 2 involving the quark 
and gluon fields in QCD, it is normally assumed that the OPE starts at dimension d = 4 (with 
the quark condensate being multiplied by the quark mass). This is supported by results from 
QCD sum rule analyses of r-lepton decay data, which show no evidence of d = 2 operators 
|15j-|17j. The unit operator in Eq.(4) has dimension d = and CqI stands for the purely per- 




Figure 3: Gluon propagator modification due to (infrared) gluons interacting with the physical 
QCD vacuum, and involving the gluon condensate. Large momentum flows only through the 
quark propagators. 

turbative contribution. The Wilson coefficients as well as the vacuum condensates depend on 
the renormalization scale. In the case of the leading d = 4 terms in Eq.(4) the /x^ dependence 
of the quark mass cancels the corresponding dependence of the quark condensate, so that this 
contribution is a renormalization group (RG) invariant. Similarly, the gluon condensate is also 
a RG invariant quantity, hence once determined in some channel these condensates can be used 



throughout. At dimension d = 6 there appears the four-quark condensate, obtained from Fig.l 
at the next to leading order (one gluon exchange) and allowing all four quark lines to interact 
with the physical vacuum (see Fig. 4). While this condensate has a residual renormalization 
scale dependence, this is so small that in practice it can be ignored. The four-quark condensate, 
while relatively small, is crucial to explain the large p(770) - ai(1260) mass splitting. In most 
applications —q'^ is chosen large enough so that the condensates of higher dimension (d > 8) 
can be safely ignored. The numerical values of the vacuum condensates cannot be calculated 




Figure 4: The four-quark condensate of dimension d = 6 in the OPE. This is responsible for the 
p — ai mass splitting. Large momentum flows only through the gluon propagator. 

analytically from first principles as this would be tantamount to solving QCD exactly. One 
exception is that of the quark condensate which enters in the Gell-Mann-Oakes-Renner relation, 
a QCD low energy theorem following from the global chiral symmetry of the QCD Lagrangian. 
Otherwise, it is possible to extract values for the leading vacuum condensates using QCD sum 
rules together with experimental data, e.g. e~^e~ annihilation into hadrons, and hadronic decays 
of the r-lepton. Alternatively, as lattice QCD improves in accuracy it should become a valuable 
source of information on these condensates. 

Turning to the hadronic sector, bound states and resonances appear in the complex energy 
(squared) plane (s-plane) as poles on the real axis, and singularities in the second Riemann 
sheet. In addition there will be multiple cuts reflecting non-resonant multi-particle production. 
All these singularities lead to a discontinuity across the positive real axis. Choosing an integra- 
tion contour as shown in Fig. 5, and given that there are no further singularities in the complex 
s-plane, Cauchy's theorem leads to the finite energy sum rule (FESR) 
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Figure 5: Integration contour in the complex s-plane. The discontinuity across the real axis 
brings in the hadronic spectral function, while integration around the circle involves the QCD 
correlator. 



f^o I 1 / 

/ ds - f{s) ImU{s)\HAD = -IT- f ds f{s)U{s)\QCD , 



(5) 



where /(s) is an arbitrary (analytic) function, Sfh is the hadronic threshold, and the finite radius 
of the circle, sq, is large enough for QCD and the OPE to be used on the circle. Physical observ- 
ables determined from FESR should not depend on sq. In practice, though, this independence 
is not exact, and there is usually a region of stability where observables are fairly independent 
of So, typically in the range sq — 1 — 3 GeV . The variation of an observable in the stability 
region is incorporated into the error of the determination. Equation (5) is the mathematical 
statement of what is usually referred to as quark-hadron duality. Since QCD is not valid in the 
time-like region (s > 0), in principle there is a possibility of problems on the circle near the real 
axis (duality violations). I shall come back to this issue later. The right hand side of this FESR 
involves the QCD correlator which is expressed in terms of the OPE as in Eq. (4) . The left hand 
side calls for the hadronic spectral function which is written as 

ImIl{s)\HAD = I'mIl{s)\poLE + I'mIl{s)\jiES + ImIl{s)\pQCDd{s - sq) , (6) 

where the ground state pole, absent in some channels, is followed by the resonances which merge 
smoothly into the hadronic continuum above some threshold sq. This continuum is expected 
to be well represented by PQCD if sq is large enough. Due to this, if one were to consider 
an integration contour in Eq.(5) extending to infinity, the cancellation between the hadronic 
continuum in the left hand side and the PQCD contribution to the right hand side, would 
render the sum rule a FESR. The performance of the contour integral in the complex s-plane is 
discussed in the next section. 



3 Finite energy sum rules 

The integration in the complex s-plane of the QCD correlator is usually carried out in two 
different ways. Fixed Order Perturbation Theory (FOPT) and Contour Improved Perturbation 
Theory (CIPT) |18j-|19j. The first method treats running quark masses and the strong cou- 
pling as fixed at a given value of sq. After integrating all logarithmic terms (ln(— s/^^)) the 
RG improvement is achieved by setting the renormalization scale to /i^ = —sq. In CIPT the 
RG improvement is performed before integration, thus eliminating logarithmic terms, and the 
running quark masses and strong coupling are integrated (numerically) around the circle. This 
requires solving numerically the RGE for the quark masses and the coupling at each point on 
the circle. The FESR Eq.(5), with /(s) = 1, in FOPT can be written as 

(-)^C2iV+2(0lO2iV+2|0)= / ds S^ - ImU{s)\HAD - S^+^ M2N+2{S0) , (7) 

Jo TT 

where the dimensionless PQCD moments M2N+2iso) are given by 

M2N+2{so) = (jv_^^ / dss'^ -ImU{s)\pQCD , (8) 

Sq jo TT 

and /?TT-n(s) is assumed dimensionless for simplicity. 
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Figure 6: Running strange quark mass in FOPT at a scale of 2 GeV as a function of sq. 
Curve (a) is for A = 330 MeV, and curve (b) for A = 420 MeV, corresponding, respectively, to 
as{M^) = 0.31 and as{M^) = 0.36. 



If the hadronic spectral function is known in some channel from experiment, then ImIl{s)\HAD = 
l7nIl{s)\DATA, and Eq.(7) can be used to determine the values of the vacuum condensates. 
Subsequently, Eq.(7) can be used in a different channel to determine the masses and couplings 
of the hadrons in that channel. It is important to mention that the correlator n(g^) is generally 
not a physical observable. However, this has no effect in FOPT as the unphysical constants in 
the correlator do not contribute to the integrals. The situation is quite different in CIPT where 
Eq.(5) cannot be used for unphysical correlators. For instance for a correlator whose physical 
counterpart is the second derivative (needed to eliminate a first degree polynomial), Cauchy's 
theorem and the resulting FESR must be written for the second derivative. In this case one has 
to use the following identity [6]- [7] 



/ dsg{s) U{s) = i ds [F{s) - F(so)] n"(s) 

Jc(\so\) Jc{\so\) 



where 



Fis) 



ds' 



ds"g{s") 



sa 



ds"g{s") 



(9) 



(10) 



and g{s) is an arbitrary function. This is easily proved by integrating by parts the right hand 
side of Eq.(9) and using Eq.(lO) to obtain the left hand side. In this case Eq.(5) becomes 



so 1 X 

ds g{s) - ImI{{s)\HAD = - 7-- 
sth TT ivrz 



C{\sq\) 



ds[F{s)-F{so)\Il"{s)\QCD 



(11) 



which is the master FESR to use in CIPT. The running quark masses and the running strong 
couphng entering n"(s) are now functions of the integration variable and are not fixed as pre- 
viously in FOPT. The running coupling obeys the RGE 

s ^^^ = /3(as) = -EI3n a^i-sf^' , (12) 

where a^ = Os/tt, and e.g. for three quark flavors /3o = 9/4, /3i = 4, /32 = 3863/384, ^^ = 
(421797/54 + 3560C(3))/256. In the complex s-plane s = so e*^ with the angle x defined in the 
interval x € (— 7r,7r). The RGE then becomes 



das(x) 

— —I 



dx 

N=0 



Y, Pn a,{xf+^ , (13) 



This RGE can be solved numerically at each point on the integration contour of Eq.(ll) using 
e.g. a modified Euler method, providing as input as{x = 0) = as(— sq). Next, the RGE for the 
quark mass is given by 



s dm{—s) 



7(a,) = - 5] 7Af a^'^' , (14) 



m ds 

M=0 

where e.g. for three quark flavors 70 = 1, 71 = 182/48, 72 = [8885/9 - 160C(3)]/64, 73 = 
[2977517/162 - 148720 C(3)/27 + 2160 C(4) - 8000C(5)/3]/256. With the aid of Eqs. (12)-(13) 
the above equation can be converted into a differential equation for m{x) and integrated, with 
the result 

m{x) = m{0) expl — i dx' N^ 7m [^3(2;')] ^ f ) (1^) 

•^° A/=o 

where the integration constant ?tt,(0) can be identified with ?tt,(so). 



4 THE LIGHT PSEUDOSCALAR CORRELATOR 



I discuss now the use of FESR to determine the values of the QCD light quark masses [6]-[7]. As 
mentioned in the Introduction, in this case the ideal current operator in Eq.(l) is the axial- vector 
current divergence. 

Mq') =^ I d^x e^"- < \T{d^A^{x)^, , 9^^4(0)1^1 > , (16) 



where d^A^{x)\\ = {rrii + mj) : 'qj{x) i'j^ qi{x) : is the divergence of the axial-vector current, and 
{i, j) are flavour indexes . To simplify the notation we shall use in the sequel rui + nij = m. The 
advantage is that the masses appear here as overall multiplicative factors, rather than as sub- 
leading power corrections like in other correlators, e.g. the vector or axial-vector correlators. 
The great disadvantage is that there is no direct experimental data beyond the pseudoscalar 
meson poles, i.e. the hadronic resonance spectral function, Im'n{s)\RES in E<l-(6) is not known 
experimentally. The only available information is that there are two radial excitations in the non- 
strange (vr) as well as in the strange (K) channel with known masses and widths. This is hardly 
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Figure 7: Running strange quark mass in CIPT at a scale of 2 GeV as a function of sq. Curve 
(a) is for A = 330 MeV, and curve (b) for A = 420 MeV, corresponding, respectively, to 
a,(M2) = 0.31 and as{M^) = 0.36. 



enough to reconstruct the full spectral function. In fact, inelasticity, non-resonant background, 
and resonance interference are impossible to guess so that a model is needed for the resonant 
spectral function. This fact, which introduces a serious systematic uncertainty, has affected 
all quark mass determinations using QCD sum rules until recently [6]-[7j. The breakthrough 
has been to introduce an integration kernel in the FESR tuned to suppress substantially the 
resonance energy region above the ground state. This kernel is of the form 



A5(s) = 1 — aQS — ai s 



(17) 



where A5(s) stands for either /(s) in Eq.(5) for FOPT, or g{s) in Eq.(9) for CIPT. The coef- 
ficients are fixed by requiring that A5(s) vanish at the peak of the two radial excitations, i.e. 
A^{M1) = A5(M|) = 0. This has the effect of reducing the resonance contribution to the FESR 
to a couple of a percent of the ground state contribution, well below the uncertainty due to a^. 
Clearly, other more elaborate choices of the integration kernel are possible. It has been found, 
though, that the simplest form above is optimal in the sense of simplicity and of achieving the 
goal of reducing considerably the systematic uncertainty from the hadronic resonance sector. 



The FESR following from Cauchy's theorem, Eq.(5), takes now the form 
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(18) 
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Figure 8: Up quark mass at 2 GeV as a function of sq for as{M^) = 0.335(0.353), or Aqcd 
365 (397) MeV, upper (lower) curve, respectively. 



where fp and Mp stand for the pseudoscalar meson pole (vr, K) parameters, and this FESR is 
suitable for FOPT. In the case of GIFT the FESR, Eq.(ll), becomes 



^/ ds 4^''''{s)[F{s)-F{so)] 

2^^ JCi\so\) 



= 2 3 Ml, A5(M|) 



1 

+ - 

vr 



ds Im ip5{s)\RES A5(s) 



(19) 



Sth 



where F{s) is defined in Eq.(lO). It should be clear that the distortion introduced by the 
integration kernel A5(s), Eq.(17), affects all three terms of the FESR as Cauchy's theorem 
remains valid as long as the kernel is an analytic function. Hence, the suppression achieved in the 
hadronic resonance contribution is compensated by corresponding changes in the pseudoscalar 
meson pole and in the QCD contributions. Since these two terms are reasonably well known, 
this is a welcome feature. 



5 RESULTS 



The light-quark pseudoscalar correlator in PQCD is known up to fifth-loop order [20], with 
the strong coupling being determined from data on r-decays |12j-|14j. as{M^) = 0.344 it 0.009 



which corresponds to a QCD scale in the MS scheme of A = 365 — 397 MeV. The handling of 
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Figure 9: Down quark mass at 2 GeV as a function of sq for as{M^] 
365 (397) MeV, upper (lower) curve, respectively. 



0.335(0.353), or Aqcd 



logarithmic quark mass singularities in this correlator requires some care, as explained in [21] - 
|22j . The leading non-perturbative contributions of dimension d = 4 are also known, together 
with higher order quark mass corrections. Complete expressions of all these contributions to 
the FESR for both FOPT and CIPT may be found in [S]-[7]. A posteriori, it is found that 
quark mass terms of order 0{m'^) and higher, as well as vacuum condensates of dimension d > 6 
are negligible on account of the integration kernel, Eq.(17). The hadronic resonance spectral 
function has been modeled by two Breit-Wigner forms normalized at threshold according to 
chiral perturbation theory [2]-[S]. The quark masses were determined with and without this 
contribution in order to gauge its impact on the final result. It turns out that allowing for a 
± 30% uncertainty in the hadronic resonance contribution impacts on the quark masses at the 
level of less than 1%. Assuming the unknown six-loop PQCD contribution to be equal to the 
known five-loop term also has an impact at the 1% level. The presence of the integration kernel, 
Eq.(17), thus enhances considerably the importance of the better known contributions to the 
FESR, i.e. PQCD and the pseudoscalar meson pole. The major sources of uncertainty are then 
the values of the strong coupling and of the radius of the integration circle sq. Results for the 
strange quark mass in FOPT are shown in Fig. 6. The width of the stability region is typical of 
FESR applications. However, in CIPT as shown in Fig. 7 this region is remarkably wide, with 
ms{2 GeV) being exceptionally stable. This is also the case for the up- and down-quark masses 
in CIPT as shown in Figs. 8 and 9. While there is no clear cut criterion to establish which 
integration technique is best, results from this application clearly favour CIPT in terms of the 
stability region. Taking into account all possible sources of uncertainty, and following a very 
conservative approach of adding them up, rather than combining them in quadrature, the final 



10 



results (in the MS scheme) are 

m„(2 GeV) = 2.9 ± 0.2 MeV , 

md{2 GeV) = 5.3 ± 0.4 MeV , 

m„, ^ !!hL±II}± = 4.1 ± 0.2 MeV , 

ms{2 GeV) = 102 ± 8 MeV . 

This is at present the most accurate determination of the light quark masses using QCD sum 
rules, with the above results in very good agreement with the most recent lattice QCD values 
[23j-[28]. With the systematic uncertainties from the hadronic resonance sector under very good 
control, a future reduction of the errors in the quark masses will mostly rely on more accurate 
determinations of the strong coupling. 
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